Frustration and Packing in Curved-Filament Assemblies: From Isometric to 

Isomorphic Bundles 
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Densely-packed bundles of biological filaments (filamentous proteins) are common and critical 
structural elements in range of biological materials. While most bundles form from intrinsically 
straight filaments, there are notable examples of protein filaments possessing a natural, or intrinsic, 
curvature, such as the helical bacterial fiagellum. We study the non-linear interplay between ther- 
modynamic preference for dense and regular inter-filament packing and the mechanical preference 
for uniform filament shape in bundles of helically-curved filaments. Geometric constraints in bun- 
dles make perfect inter-filament (constant spacing, or isometric) packing incompatible with perfect 
intra-filament (constant shape, or isomorphic) packing. As a consequence, we predict that bundle 
packing exhibits a strong sensitivity to bundle size, evolving from the isometric packing at small 
radii to an isomorphic packing at large radii. The nature of the transition between these extremal 
states depends on thermodynamic costs of packing distortion, with packing in elastically-constrained 
bundles evolving smoothly with size, while packing in osmotically-compressed bundles may exhibit a 
singular transition from the isometric packing at a finite bundle radius. We consider the equilibrium 
assembly of bundles in a saturated solution of filaments and show that mechanical cost of isomorphic 
packing leads to self-limited equilibrium bundle diameters, whose size and range of thermodynamic 
stability depend both on condensation mechanism, as well as the helical geometry of filaments. 

PACS numbers: 



I. INTRODUCTION 

Molecular filaments and their assemblies underlie the 
function of a diverse range of materials, from intra- and 
inter-cellular matter to high-strength synthetic materi- 
als. Condensed, parallel arrays, or bundles, of cytoske- 
latal filaments like f-actin and microtubules provide both 
mechanical reinforcement to cells and play crucial roles 
in vital cellular processes like adhesion, division and lo- 
comotion P, Rope-like assemblies of extra-cellular 
protein filaments — for example, collagen, cellulose or 
fibrin — are also primary structural elements in-between 
cells, in plant and animal tissue [!, H[. As filaments are 
densely-packed and hence strongly interacting in bundles, 
the structure of intermolecular organization underlies the 
collective dynamic and mechanical properties of these as- 
semblies, critical for their robust function in living sys- 
tems. In the simplest models of this material prototype, 
the ground state of bundles of straight filaments is uni- 
versally a hexagonally-packed array of parallel filaments, 
which maximizes both the number of inter-filament con- 
tacts as well as the density within the bundle. As a con- 
sequence, filament bundles forming under the influence of 
cohesive inter-filament forces @ as well as osmotic pres- 
sure (depletion forces) @ are observed to form hexagonal 
cross-sectional order. Notably both the (straight) fila- 
ment shape and the inter-filament spacing are uniform 
throughout a hexagonal bundle, which implies that both 
length and lateral size of self-assembled bundles is ther- 
modynamically unlimited 0, H| ■ 

In this article, we consider the more complex structure 
and thermodynamics of bundles formed by intrinsically- 



curved filaments, a distinction which we show frus- 
trates homogeneous inter- and intra-filament geometry 
of straight, hexagonal bundles. For many common biofil- 
amcnts — f-actin or microtubules — the filament back- 
bone is straight in the absence of external stress or ther- 
mal fluctuations. There are, however, important excep- 
tions, where the backbone filament geometry has a pre- 
ferred, or natural, curvature, deriving from the shape 
and interactions between globular subunits Q. Most 
widely known among this class of intrinsically-curved fil- 
aments is the bacterial flagella [10|, a corkscrew shaped 
protein filament — helical radius r = 0.2/im; and pitch 
p = 2.2/im; diameter s» 20 nm, shown in Fig. [T] (b) — 
that propels common microorganisms like E. coli and 
salmonella. Notably, bacteria swim, not by the actua- 
tion of single filaments, but by the screw-like rotation of 
coherent bundles of multiple flagella [ll|, held together 
through a combination of contact and hydrodynamic in- 
teractions. More recently, new classes of curved-filament 
assemblies have been discovered as components of bacte- 
rial "cytoskeleton" Q. One such filament MReb, consid- 
ered to be a prokaryotic analogue to f-actin, is observed to 
associate with inner-cell wall of bacteria in bundles that 
wind helically around the cell body, potentially shaping 
the growth of rod- like bacteria [12|, Qij- A second fila- 
ment, crescentin, also forms bundles in vivo fl4| whose 
helical structure is believed to drive the curved and heli- 
cal shapes of bacteria like Caulabacter crescentus [15j. 

Unlike the case of straight-filament assembly, bundles 
of curved-filaments are subject to an unavoidable geo- 
metric conflict between inter-filament packing and the 
natural curvature of filaments. Simply put, it is not pos- 
sible to evenly space nearest-neighbor filaments through- 
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FIG. 1: Preferred (stress free) geometry of helical filaments 
shown in (a), where tan 9h = 2-wr/p. In (b), florescence- 
microscopy images of a single bacterial flagella (inset) and 
a bundle of (~ 10) flagella formed in a neutral polymer sus- 
pension, scale bars are 5 fim (courtesy of Z. Dogic, Brandeis 
University). The two extreme cases of bundle packing of he- 
lical filaments: (c) the isometric (constant spacing) packing; 
and (d) the isomorphic (constant shape) packing. The orange 
filament on the surface of the bundles highlights the subtle 
changes of filament shape required by the isometric packing 
in (c), while the lateral spacing in (d) is such that filament 
diameters do not overlap in the bundle interior. 



out the ordered cross-section of the bundle while main- 
taining uniform shape of all filaments. This conflict is 
illustrated by the bundles of helical filaments in Fig. Q] 
(c) and (d), which we denote, respectively, as isomet- 
ric (constant spacing) and isomorphic (constant shape) 
packings. 

The perfect inter-filament packing — constant center- 
to-center filament distance — in the isometric bundle 
constrains the relative geometry between filaments so 
that the inter-filament separation, A, between any fil- 
ament pair is relatively parallel to the filament tan- 
gents [16J. This means that in isometric bundles all fila- 
ment backbones are normal to the common plane of lat- 
tice packing, and A does not twist around the backbones 
along the length of the bundle [Fig. Q] (c)]. Importantly, 
this packing constrains the relationship between the cur- 
vature k and torsion r of a "reference" filament and the 
shape of a filament at separation A, whose respective 



where N is the normal to reference filament. Eq. ([I]) 
is the filamentary analogue of the "curvature focussing" 
properties of evenly-spaced, layered materials (i.e. smec- 
tics) (l7l - i20j . embodying the fundamental frustration of 
the lateral assembly of filaments in densely-packed bun- 
dles. On one hand, driving forces for dense-filament 
packing or cohesive filament generically penalize devia- 
tions from perfect lateral packing in a bundle. However, 
when these filaments are intrinsically curved, meaning 
that k =^ in their stress free configurations, perfect lat- 
eral packing requires the filament geometry to vary along 
normal directions (N • A ^ 0) in the bundle cross sec- 
tion, incurring a mechanical cost for change of filament 
shape. Thus, in bundles of intrinsically-curved filaments, 
isometric packing is not compatible with constant fila- 
ment shape, a property achieved in the isomorphic pack- 
ing of Fig. [T] (d). In the isomorphic state, as we show 
below, the non-local nature of filament contact requires 
inter-filament spacing to vary along the bundle, straining 
either interactions that prefer constant spacing or ther- 
modynamic forces favoring maximal density. The mutual 
frustration between isometric and isomorphic packing im- 
plies that the state of inter-filament packing in curved- 
filament bundle is determined not only by the lateral 
forces that drive filaments in to tightly-packed assem- 
blies, but also by the mechanical costs of filament dis- 
tortion. Here, we show that interplay between geometry 
and mechanics in curved-filament bundles has a profound 
influence on the structure and dimensions of assembly, 
giving rise to thermodynamic behavior not exhibited by 
assemblies of unfrustrated, straight filaments. 

In this paper, we study the frustration between con- 
stant spacing and constant shape in bundles of helically- 
shaped bundles, whose backbones have preferred values 
of curvature and torsion, kq and tq . Helical filaments are 
prototypes of the frustration in curved-filament assem- 
blies, whose degree of frustration can be tuned continu- 
ously from zero for straight filaments (kq/tq <C 1) and 
the maximal frustration of ring filaments (ko/tq 3> 1). 
In ref. Q, it was shown that the mechanical costs of 
filament shape change in isometric, helical filament bun- 
dles grow rapidly with lateral size of bundles, ultimately 
limiting the lateral radius of bundles that are unable to 
relax the constant-spacing constraints. In the present 
study, we develop a continuum theory to analyze the 
non-linear interplay between the distinct costs of inter- 
filament and intra-filament distortions in helical-filament 
bundles which relax both shape as well as the lateral 
packing. We study two models for the thermodynamic 
cost of inter- filament distortion: 1) an elastic model for 
disruption of inter-filament cohesion; and 2) the osmotic 
excluded- volume cost of bundles formed in concentrated 
solutions of osmolytes excluded from the filament pack- 
ing. Analyzing a class of low-energy bundles that span 
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continuously from the limiting isometric to isomorphic 
bundles, we show that bundle packing evolves from per- 
fect inter-filament spacing in the limit of vanishing bun- 
dle diameter to constant hiament shape in the limit of 
large diameter. While inter-filament packing in cohesive 
bundles is distorted for any finite radius, we find in the 
limit of weak resistance to filament torsion, osmotically- 
condensed bundles remain isometric over finite range of 
small bundle size. Finally, the consider the equilibrium 
size of bundles in the canonical ensemble, showing that 
the mechanical cost of shape change in small bundles 
stabilizes finite size bundles when the surface energy of 
bundles is below a threshold value, which is shown to 
increase with kq/tq for both bundling mechanisms. 

The remainder of this article is organized as follows. 
In Sec. II. A we introduce the geometry of helical fila- 
ment bundles, as well as the mechanical cost for filament 
distortion and the thermodynamic costs for distortion of 
the packing in Sec. II. B. In Sec. III. A we analyze the 
variation of bundle packing with lateral size of bundles 
in both cohesive and osmotically-driven bundles, and in 
Sec. III.B we consider the equilibrium radius of bundles 
in solutions of freely-associating filaments above satura- 
tion. Finally we conclude in Sec. IV with a discussion 
of the implications on our study for bundles of flagellar 
filaments. 



II. STRUCTURE AND MECHANICS OF 
INTER-FILAMENT PACKING IN BUNDLES 

We consider two models of helical-filament bundles, 
which differ in terms of the driving forces for assembly. 
In the first model, assembly is driven by short-range, 
adhesive interactions between neighbor filaments in the 
bundle, while the second model considers densely packed 
bundles to form under the influence of osmotic pressure, 
that works against the total volume of the bundle inte- 
rior. In both models, an increase of inter-filament pack- 
ing forces (adhesion or pressure) drives the bundles away 
from the isomorphic state preferred by unbundled fila- 
ments towards the isometric packing. To describe this 
transition, we consider a class of bundle configurations 
that may be continuously transformed from the isomor- 
phic into the isometric state and vice versa. 



A. Geometry of filament packing 

We consider bundle configurations described by a pack- 
ing the cross-section around a central filament, whose 
center is described by, Ro(s). The backbone of the cen- 
tral filament is described by a tangent, normal or bi- 
nomial, the Frenet frame {T, N, B}, curvature k and 
torsion r, deriving from the relations, 

<9 s Ro = T; d s T = kN; <9 s N = -rB. (2) 



The curvature and torsion of this filament define the he- 
lical angle 9h, 

tan Oh — k/t, (3) 

shown in Fig. [T] (a). In the bundle cross-section, fil- 
aments are labeled by coordinates (x u ,x v ). The vector 
x = x u u + x v v maps the position of the central filament 
to the filament at {x u ,x v ) in the plane spanned by the 
basis vectors, u and v, 

R(x u ,x v ,s) = Ro(s) +x. (4) 

We refer to the plane normal to p = u x v as the packing 
plane, and the orthonormal frame {u, v, p} as the packing 
frame. We assume a regular, hexagonal lattice packing, 
i.e. x u — a(n + m/2) and x v — amyZ/2 where m and 
n are integers labeling occupied lattice positions and a 
lattice spacing. 

Due to the constraints of packing, the geometry of fila- 
ments varies through the cross-section of a bundle. Given 
our construction, we describe the variation in terms of the 
functions T(x), /c(x) and r(x) characterizing respective 
values of backbone orientation, curvature and torsion for 
filaments at x. We focus on the limit of narrow bun- 
dle cross-sections, such that k|x| <C 1 and r|x| <C 1. 
The shape of the filament at x derives from both the 
shape of Ro(s) as well as the tilt and rotation of the 
frame {u, v, p} with respect to the geometry of the cen- 
tral curve at x = 0. The rotation of filament positions in 
the packing plane along the arc length of the bundle is 
described by the transformation, 

u = cos(51s)u' + sin(f2s)v' 

v = - sin(Os)u' + cos(Qs)v', (5) 

where 

u'=N; v' = sin6»T + cos6>B; p = cos 9T - sin 0B . (6) 

Here, 9 is the tilt angle of the packing plane, the pack- 
ing angle, with respect to the central tangent, i.e. cos 9 = 
T-p. Thus, SI and 9 constitute two independent parame- 
ters specifying the geometry of the packing. We consider 
bundles that are homogeneous along their length so that 
both SI and 9 are constant. 

The relationship between geometry (i.e. orientation, 
curvature and torsion) of the central filament to the ge- 
ometry of a filament at x derives from rotation of the 
{u,v,p} frame along the bundle, which follows straight- 
forwardly from eqs. (|5|) and © 

d s u = (SI + Oq)v — Sl p cos(£ls)p 

d s v = — (SI + Qq)u + Sl p sin(f2s)p 

d s p = Sl p u' = QjJcos(f2s)u + sin(fis)v], (7) 

where 

SIq = — v' ■ s \i' = — k sin 9 -\- t cos 9 

Sl p = — p • B s ii f = k cos 9 -\- t sin 9. (8) 
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From eq. we note that the net rotation of fila- 
ment positions around the center of the cross section, 
v • d s u = (SI + Slo) is the simply the sum rotation of u' 
and v' around p, Slo, and the rotation of u and v rela- 
tive the to that coordinate system, SI. In our analysis, we 
focus on the set of low-energy states that span the two- 
limiting packings, isometric and isomorphic. We show in 
the appendix that vanishing torsion of filament positions 
in the packing plane (i.e. SI + Slo = 0) is a necessary 
condition for both constant filament shape and constant 
filament spacing, and we therefore, expect this condi- 
tion to also describe the low-energy states intermediate 
to these extreme cases. That is, we restrict our analysis 
to 

n = -O , (9) 

so that there is no rotation of the filament position in the 
packing plane along the bundle. 

The variation of filament geometry throughout the 
packing plane derives from a calculation of the Frenet 
geometry based on the construction of eqs. ©-(|H|) (sec 
Appendix). For backbone orientation (tangents) we find, 
first order (in k|x|) correction T(x) — T(0) ~ ST± 

8T ± = <9 s x — T(<9 s x ■ T) 

= sin6%(N-x)B, (10) 

which is manifestly perpendicular to T. The variation of 
backbone curvature and torsion are described by 

k(x) ~ k + Sk (kN ■ x); r(x) ~ r + St (kN • x) (11) 

where 

Sk = Sl p (2 cos 9 — Sl p /k) 

= k(1 - sin 2 6>/sm 2 6> p ) (12) 

and 

St = Sip ( cos # t/k + sin 9 — sin^Sl 2 /*; 2 ) 
= kcos 2 9(ta,n9 cot 9 P — 1) [tan0 — cot9 p 

- cos 9 sin 9 (cot 9 p + tan 9) 2 ], (13) 

where the common root of expressions occurs for packing 
angle, 

9 P = - arctan(K/r) = -6 h . (14) 

Notably, shape change of filaments (to lowest order) is 
distributed anisotropically in the packing plane: curva- 
ture and torsion are maximally distorted for filaments 
along x || N, while no shape change is required for fila- 
ments at x 1 N. 

This construction encompasses two extremal limits of 
the packing geometry: 

Isomorphic Packing: Sl p = - Inspection of eqs. (|10[) 
- (113[) immediately reveals these conditions achieve uni- 
form orientation and shape the packing plane, ST — Sk = 



St = 0. Solution of eq. (|8]) gives the inclination and ro- 
tation rate, SI, the isomorphic packing, 

9 = 9 p ; n = -Vk 2 +t 2 . (15) 

In this packing, the packing plane is fixed (say, the x — y 
plane in the lab frame) and remains normal to the pitch 
axis of the helices (z direction). Though filaments are 
parallel in the packing plane, inter-filament distances in 
this plane are not the true distances between the helical 
curves. That is, u and v are not relatively parallel to the 
backbone, T. Below we show that this implies variation 
of inter-filament spacing the pitch axis of the helices in 
isomorphic bundles. 

Isometric Packing: 9 = - In the isometric packing, 
filaments are parallel in the packing packing plane, p = T 
and there is no twist of the packing around T. Because u 
and v are relatively parallel to T, distances in the pack- 
ing plane are true inter-filament distances, and hexagonal 
packing in this plane implies uniform neighbor spacing 
along the length and throughout the cross section of the 
bundle. While spacing is constant, shape varies through- 
out the packing plane according to 

Sk = k; St = r, (16) 

which derive from eqs. (jT2"j) and (TT5]) for Oq = ~ T an d 
ftp — k, consistent with the shape variation of eq. ([T]). 

B. Mechanics of filament distortion 

In this section, we describe the mechanical cost of dis- 
tortions of filaments away from a preferred helical shape 
in bundles, characterized by constant values of backbone 
curvature kq and torsion r . We adopt the Kirchoff-Love 
theory of filament mechanics, where distortions are cap- 
tured the rotation of a material frame along the back- 
bone [2 it] . For each filament, labeled by position x, 
this frame is described by the triad {ei(x), e2(x), e3(x)} 
where e3(x) = T(x) is along the backbone tangent and 

ei(x) = cos 0(x)N(x) + sin0(x)B(x) 

e 2 (x) = - sin 0(x)N(x) + cos 0(x)B(x), (17) 

describe the position of the filament cross section, and <j) 
represents an internal, rotational degree of freedom of the 
filament. The geometry of filament backbone is locally 
described by three rotation rates, 

de 3 de 3 dei . . 

Ki=ei - ——; k 2 = e 2 • -7-; to = e 2 • -3—, (18) 
as as as 

where for simplicity we have suppressed the x dependence 
of each quantity. The first two rates, k\ — ncoscj) and 
K2 = —Ksincf) describe the bending of the filament back- 
bone projected onto the two material directions in the 
filament cross section, while the last rate describes the 
rotation, or torsion, of cross-section around the filament 
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FIG. 2: On the left, three examples of the helical filament packings considered in our model, spanning from isometric (6 — 0) 
to isomorphic (6* = Q p ). The orientation of the "packing plane" of hexagonal order is shown in each. In each bundle, orange 
filaments highlight a single row of filaments, whose bending and inter-filament packing geometry is depicted on the right for the 
same three packing angles. On the top right, the color variation highlights local deviations from preferred curvature, |re(x) — k| 2 , 
ranging from blue for undistorted to yellow for highly-distorted curvature. In the bottom right, the color variation highlights 
strain of interfilament distance along the row, deriving from the local tilt in the packing plane along the neighbor direction 
according to eq. (|21[) . where blue indicates no strain and yellow indicates high strain. 



backbone, uj = t + d s cj>. In this framework, a helical fila- 
ment is described by a preferred constant value of twist 
To, curvature kq, and a preferred direction of bending, 
which we choose to be ei. 

For simplicity, we focus our attention on the limit of 
"easy twist" filament mechanics [8|, which can be justi- 
fied strictly for the case where the bending modulus of 
the filament B is much larger than the twist modulus 
of the filament C. In this limit, for any modulations of 
«(x) and t(x) the torsional state of the filament adjusts 
to maintain alignment of the normal to ei, the direc- 
tion of preferred curvature. Explicitly, 4> — d s <f) = and 
the mechanical cost of shape change for a filament at x 
becomes simply, 

£ fil (x) = i J d S (x){s[ K (x)- Ko ] 2 + C[r(x)-ro] 2 }. 

" , ^ (19) 

To describe a filament bundle, we consider uniform areal 

density of filaments p Q and integrate -Efli(x) over circular 

bundle cross section of radius R. Using cq. (fTTj) we find 

the total mechanical energy for a bundle N filaments 

£mech B 2 C. , 2 

+ ^Sk 2 (kR) 2 + ^5t 2 {kR) 2 , (20) 
8 8 

where L is the filament contour length. In deriving eq. 
(|20[). we have assumed that shape of the central filament 
remains close the preferred shape for sufficiently narrow 
bundles, specifically, \n — Kq\R ~ \t — tq\R ~ (fto-R) 3 , 
which follows from minimization of eq. (|20p with respect 
to K and r. In the remainder of the article, we drop the 
distinction drop the distinct between the preferred helical 
shape and shape of the central filament, and take k = kq 
and t — tq. In Fig. [5] we show graphically the spatial 



variation of bending energy density in bundles ranging 
from bundles ranging isometric to isomorphic packing. 



C. Thermodynamics of inter- filament packing 

We consider two models of the thermodynamic cost of 
inhomogeneous filament packing: 1) a continuum elastic 
model of cohesive filament bundles; and 2) the osmotic 
cost of excluded volume in densely packed bundles. In 
both models, free energy costs arise as the consequence of 
the non-local nature of inter-filament contact. Consider 
two filaments with backbones centered at x and x + <5x 
in the packing plane. Assuming that ko < 1, filament 
backbones may be approximated locally as straight seg- 
ments of orientation, T(x) ~ T(x+i5x) for nearby points. 
Rather than the in-plane separation (foe, interactions be- 
tween homogeneous filaments are natural functions of 



<frc-T(5x-T), 



(21) 



which is the distance of closest approach between the 
filament centers H3- Since |A ± | 2 = |5x| 2 - (<Jx • T) 2 
when filaments are locally tilted into the packing plane, 
inter-filament spacing is necessarily reduced relative to 
the packing plane distance. From eq. (|10p . the in-plane 
components are 

T„(x) = sin(Qs)sin0[l + cos0Qp(N-x)] 

T„(x) = cos(fis)sin0[l + cos0fi p (N-x)], (22) 

which shows that for 8 ^ in-plane tilt (in the packing 
frame) varies along bundle, as well as within the bundle 
cross-section (for £l p 7^ 0). The interfilament strain gen- 
erated by tilt in the packing plane is depicted in Fig. [5] 
graphically for bundles ranging isometric to isomorphic. 
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Below we find that inhomogeneous patterns of tilt, lead 
directly to thermodynamic costs associated with distor- 
tions away from the isometric state. 



1. Elastic cost of inter-filament distortions 

In this section we derive a model for filaments con- 
densed by the finite-range attractive interactions, where 
the energetic cost of small changes in local inter-filament 
spacing |Aj_| are described by a non-linear continuum 
elastic theory of columnar materials [HJ • The role of non- 
linear coupling between strains and tilts has been previ- 
ously discussed and analyzed extensively in the context of 
inhomogeneous packing in chiral (twisted) filament bun- 
dles @, H, HI, HH • The elastic energy has the form, 



Selast = ^ J dV[\{u kk f + 2(1 UtjUij] . 



(23) 



where A and p, are the Lame constants corresponding to 
the isotropic elastic response of a hexagonally-ordered 



material 24|. Here, the 2-component strain, de- 
scribes changes in inter-filament spacing of nearby fila- 



ments, |Aj 



relative to initial separation |Al| 2 



5(\Aj 



(24) 



where <5x is the two-component in-plane separation be- 
tween filaments in the reference state. We consider as 
a reference state, filaments initially oriented normally to 
the (u,v) packing plane, so that indices i and j sum over 
u and v directions and |Aj_| 2 = |<5x| 2 . The form of the 
elastic strain follows, 



5(* 



djUi 



(25) 



which is sensitive only to changes in inter-filament spac- 
ing perpendicular to backbone orientation T [3l| . Im- 
portantly, eq. (1251) shows that inter-filament spacing, 
and consequently, inter-filament cohesive energy changes 
both with shifts of relative in-plane filament positions, as 
described by 



djUi), 



(26) 



and with deformations that tilt backbones into the pack- 
ing plane, as described by — TjTj/2. Note that while 
these latter terms vary along the bundles according to 
eq. (|22[) , in the packing frame the relative positions of 
filaments are fixed along s, according to eq. (j4]). This 
implies that the net effect of in-plane tilt, when aver- 
aged over the length of the bundle, is to generate an 
isotropic compression, which in turn, generates a tensile 
response in the packing. To offset the costs of tilt-induced 
compression, we consider an isotropic, in-plane strain by 
u% = (a + f3r 2 )8ij where a and /3 are coefficients param- 
eterizing the isotropic dilation of the packing at a radius 



r = |x| from the central filament. Inserting into eq. (f23|) 
and minimizing over a and /? we find, 



E r 



last 



V 



M • 4 
— sin 



8 + 0[(nRY 



(27) 



where again we retain terms to lowest order in kR and 
PoV/L = N where po is the in-plane density of filaments. 



2. Osmotic cost of excluded volume 

In this section, we derive the cost associated with vol- 
ume change of filaments bundled under the influence of 
an osmotic pressure, II. This model is relevant to the case 
of filaments in a solution of high-molecular weight poly- 
mers, whose size excludes them from the interstitial space 
between densely-packed filaments in condensed bundles. 
The osmotic free-energy cost of a bundle is simply, 



Fry 



— n v m 



(28) 



where V cyi is volume within of the bundle. Here, we con- 
sider a "hard tube" model of filament packing, which re- 
quires that the local distance of closest approach between 
any filament and its neighbor in the packing satisfies, 



I A, I >d, 



(29) 



where d is the filament diameter. As described above, 
local tilt of filaments into the packing plane frustrates 
the uniform spacing, hexagonal close-packing, requiring 
a distortion of the bundle cross-section. To avoid overlap 
along nearest neighbor direction a (unit vector in pack- 
ing plane) we assume that the lattice spacing, a, in the 
packing plane expands to maintain, 



! (1 



T(x)] 5 



(30) 



Since the bond direction a is fixed in the packing-plane 
while in-plane components of T(x) rotate along the bun- 
dle according to eq. (I2"2"j) , this condition must hold at x 
for all neighbor directions a, including the in-plane tilt 
direction v'. Neglecting the O(nr) correction to T(x) 
this gives us the packing-plane lattice dimension, 



dsec9, 



(31) 



with the areal density of filaments in the packing plane 
of p{0) — pocos 2 8, where po — 4/\/3d~ 2 . The cross- 
sectional area of the packing plane is simply, 



-4 P ack = Np 1 sec 2 



(32) 



The total volume excluded by the core of bundle is equal 
to the volume swept out by translating the cross-sectional 
area of the packing plane along the contour of the central 
filament Rq. Since the packing plane is inclined by an 
angle 6 with respect to the tangent of Ro , this excluded 
volume is simply, 



V = A 



pack 



L cos 6, 



(33) 
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and the osmotic cost forming the bundle becomes 



isomorphic 



NL 



Upo sec 9. 



(34) 



Notably, the volume-expansion cost for osmotically con- 
densed filaments grows more quickly at small 9 than for 
clastic distortion model. 



III. OPTIMAL PACKING AND SIZE OF 
HELICAL FILAMENT BUNDLES 

We now analyze the structural transition of helical fil- 
ament bundles of fixed radius in terms of the two models 
that penalize distortions of inter-filament packing, de- 
scribed by -Eciast and F osm derived in the previous section. 
We determine the equilibrium variation of packing angle 
9, which serves as an order parameter, with bundle size, 
filament shape and respective costs of inter- and intra- 
filament distortion. Specifically, we show that in the limit 
of strong packing forces, bundles approach the isometric 
state 9 — > as radius vanishes, while in the opposite 
limit of large mechanical costs for intra-filament distor- 
tion, bundles approach the isomorphic state, 9 — »• 9 p . 

First, we analyze the state of packing in both models 
for a given R, k and r, and then consider the equilibrium 
size of freely-associating filament-bundles in solution. 




R e \/ sin 8 P 



FIG. 3: In (a), the variation of optimal packing angle 
for cohesive bundles with radius of the "easy twist" limit 
(/? — > 0). The finite torsional stiffness case (/3 = 1) is 
shown in (b), for a sequence of packing angles: sin# p = 
0.125,0.25,0.375,0.5,0.625,0.75,0.875,1. 



A. Cohesive bundles: elastically-constrained 
packings 

The free energy cost of filament distortion and inter- 
filament forces in cohesive bundles is described by 



Rcoh — F mccn -f" E, 



clast , 



(35) 



where the forms of these functions are given be eqs. (|20l 
and (|27p. respectively. Defining the dimensionless mea- 
sures of curvature and torsion variation, 



5k — k 1 Sk: St = k St, 



(36) 



and grouping the ^-dependent terms into f co h(9) (free 
energy per filament) 



fcoh{9) _ Bk*R 2 



2 „:„2 , 
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{|^(6»)| 2 +/?cot 2 6'p| ( 5f(6l)| 2 +i? ol -sin 



we note the state of packing (as described by 9) is deter- 
mined by following three dimensionless parameters for 
cohesive bundles 



tan# p = -k/t; (5 = C/B; R cl = R/R cl [R cl = ^{jf^ ' 



2k 2 

_ (38) 

The first of these, tan# p , is a simply measure of central 
filament shape, while the second, /3, is simply an intrin- 
sic property of the mechanics of the filaments. The last 



of these parameters, -R e i, is a dimensionless measure of 
bundle size that characterizes the relative costs inter- vs. 
intra-filament distortions. Hence, this parameter plays a 
role similar to the Foppl-von Karman number for prob- 
lems involving confined elastic membranes. For small 
bundles, R c \ <C 1, the inter-filament elastic energies are 
large in comparison to the cost of filament shape change, 
while the inverse is true for large bundles, where .R c i 3> 1. 

To understand the size-dependence of packing we con- 
sider two cases for (3. First, we analyze the extreme "easy 
twist" limit, where f3 — > 0, and the minimization of / co h 
with respect to sin 9 is analytically tractable, resulting in 

sin 2 6Ljt 2 , \i/2 



V sin 6>„ + R z , / 



R e \/ sm9 p for 7? e i < sin 2 L , 
sin 6 p for i? cJ > sin 2 9p ' 

(37) 

As shown in Fig. [3] (a) this behavior demonstrates that 
cohesive bundles adopt the extremal isometric or isomor- 
phic packings only in the R — ¥ (single filament) and 
R — > oo (bulk assembly) limits where the respective costs 
of inter- and intra-filament distortion dominate. 

Extending the analysis to the case of a finite twist elas- 
tic cost, j3 0, we consider the small s'm9 limit of eq. 
(|37jl. finding optimal packing for narrow bundles, 



lim sin<2 coh (/3) ~ (R 2 cl p cot 3 9 p /2) 



1/3 



(40) 



8 



which indicates a much more rapid change in packing 
angle with size for small bundles, ~ i? 2 / 3 , due to finite 
torsional stiffness. The change in power-law dependence 
derives from the fact that in the small angle limit the 

5k 2 (9) - <5k 2 (0) sin 2 9, whereas 5t 2 {6) - St 2 (0) ~ 

— sin 9. Hence, the "torque" acting on packing angle gen- 
erated by bending vanishes in the isometric state (9 = 0), 
while the torque generated by torsional cost remains fi- 
nite for finite /3 cot 3 9 p . The singular, power-law evolu- 
tion of sin# con (/3) at small bundle size is also strongly 
dependent on filament geometry. As filament geometry 
approaches shape, 9 V — > 0, the small-i? evolution from 
isometric to isometric bundles becomes infinitely sharp, 
indicative of the negligible disruption of inter-filament 
packing in this limit. In the opposite limit of large tilt 
angles 9 p —> tt/2 the cost of filament torsion in the iso- 
metric state becomes negligible, and the dependence of 
packing angle on bundle size (for narrow bundles) ap- 
proaches the linear evolution of cq. (HE))) . The numerical 
solution for 9 co h{(3 = 1) is shown in Fig.[3](b) for several 
filament shapes, ranging from the rapid jump between 
isometric to isomorphic bundles for nearly straight fila- 
ments (k/t <C 1) to the more gradual evolution of pack- 
ing in highly curved filaments (k/t ^> 1). 

B. Osmotic bundles: pressure-constrained packings 




FIG. 4: In (a), the variation of optimal packing angle for 
osmotically-condensed bundles with radius of the "easy twist" 
limit (P — » 0), shown here for sm9 v = 0.5. The finite torsional 
stiffness case (j3 = I) is shown in (b), for a sequence of packing 
angles: sin 6 P = 0.125, 0.25, 0.375, 0.5, 0.625, 0.75, 0.875, 1. 



The fre- energy cost of distortion of filament shape and 
inter-filament forces in osmotically-condensed bundles is 
described by 



F n 



(41) 



where the forms of these functions are given be eqs. (|2H|) 
and respectively. Following a similar analysis to 

the above for cohesive bundles we group the ^-dependent 
terms into / osm (0) 



f osm (9) Bk 4 R 2 



{\5R(9)f 



-13 cot 2 6 p \Sf(6)\ 2 +AR- 



For pressure-induced bundles, the packing is determined 
by following three dimensionless parameters 



tanflp = -k/t; (3 = C/B: 

R os = R/R os (i? os = y/2U/ (Bk' 



Po 



(43) 



The first of these two parameters, reflecting intrinsic fil- 
ament properties, are appear identically to the analysis 
of cohesive bundles eq. (|3"8")) . where as the dimensionless 
measure of bundle size, -R s, captures the relative cost of 
bundle volume change to filament distortion, including 
the dependence on osmotic pressure. 

Again, we first analyze equilibrium packing in the j3 —¥ 
limit. The minimization of / OS m is further simplified by 
noting sin 2 9 < 1 which allows us to approximate sec 9 ~ 



1 + sin 2 9/2 + 3 sin 4 9/8. Minimizing with respect to 9 in 
this limit we find, 



sin 



,09 -> 0) = 





sin 2 S p (Rl 



i? 2 a +3sin 4 6 p /2 



for R 
for R n 



< sin 9 P 
> sin 9 p 

which is shown in Fig. 2] (a). Thus, in comparison to 
cohesive bundles, we find a more abrupt change a packing 
with bundle size for osmotic bundles. Notably we find a 
second-order transition from the isometric state (9 = 0) 
) occurring at a critical size R os = sin9 p . Near to this 
jcrhical value, for R os > sin^p, we find 9 osm ~ \R OB — 
Wrwp] 1 / 2 indicating a high sensitivity of packing near the 
critical size. As in the case of elastic restoring forces, 
we find that large bundles R os 3> sin 9 asymptotically 
approach the isomorphic packing 9 osnl — > 9 p . 

Following the analysis of finite torsional stiffness ((3 ^ 
0) above for cohesive bundles, we analyze the small sin 9 
limit of eq. (|4"!2)) to determine the small diameter evolu- 
tion of packing, 



_lim sin0 osm (/3)~i? 2 s /3cot 3 p /2. 



(45) 



Hence, for finite torsional stiffness, bundles are distorted 
from the isometric packing for any non-zero bundle ra- 
dius, deriving again from the non-zero torque exerted on 
the isometric bundle by torsional mechanics of filaments. 
Fig. 0] (b) shows the evolution of packing angle with bun- 
dle radius for /3 = 1. Notably, as we found for cohesive 
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bundles, when the helical angle approaches 9 p —> n/2, 
effects of torsional stiffness at small packing angles van- 
ish, and we find that bundles remain isometric for a finite 
range of small bundles, exhibiting a second-order transi- 
tion to 9 os at a critical size .R s- 



C. Self-limited bundle size 

Having analyzed the size- variation of optimal packing, 
we now consider the equilibrium size of bundles formed 
in concentrated filament solutions. In a previous study of 
helical filament bundles 0, which assumed the packing 
remains locked to the isometric limit, it was shown that 
the thermodynamically optimal size of bundles remains 
finite over a broad range of cohesive energy due to the 
prohibitive cost of filament shape change in large bun- 
dles. Presently, we generalize this analysis by considering 
the relaxation of filament distortion through adjustment 
of bundle from the constant-spacing towards constant- 
shape packing. 

We consider solutions in the canonical ensemble, where 
the filament number and volume are constant. We fur- 
ther assume that all but a negligible number of filaments 
belong to bundles (i.e. far above the saturation point) 
and consider a size distribution that is sharply peaked 
around radius R. In this case, an additional-free energy 
cost enters, associated with the surface energy, £, and 
the exposed bundle surface area , 



(46) 



where we assume L> J? such that the contribution from 
the bundle ends is negligible. 

To analyze the evolution of bundle size distribution 
with S, we assume |<5/ci?| and \StR\ remains sufficiently 
small that the central filament locks into the preferred 
filament geometry, k = kq and tq — r. For the purposes 
of simplifying the analysis we focus on to the limit of easy 
twist where /3 = 0, though it is straightforward to extend 
the more general analysis of f3 ~ 1. The free energy per 
filament of bundles of size R may be written as, 



Ftot(R) 2Ep 1 -Fbuik(-R) 



NL 



R 



NL 



(47) 



where Ftmik represents the bulk costs associated with 
packing distortion (optimized over packing angle, 9). 

In both models, the cost associated with the com- 
petition between inter- and intra-filament packing are 
monotonically increasing functions of R (see Fig. [5]). 
For small radius — R <C R c \ sin 2 9 P for cohesive bun- 
dles and R < R os sin 9 P for osmotic bundles — bundles 
are packed isometrically and the mechanical cost asso- 
ciated with filament distortion grows as Fb u ik/(NL) ~ 
k a R 2 . In the opposite large bundle limit, bulk cost 
of packing in both bundles (per filament) saturates as 
the bundle approaches the isomorphic state 9 — > 9 pi 
where F^J(NL) = fip^ 1 sin 4 6 V and F^/{NL) = 



1.4 

4 i.o 

£.0.6 
^ 0.4 
0.2 




R/(R os smO p ) 



FIG. 5: In (a) and (b), the respective free energy of cohe- 
sive and osmotically-condensed bundles of helical filaments 
(sin 9 P — 1 /2) as function of bundle radius for a range of sur- 
face energies as indicated in the plots. Both condensation 
mechanisms show finite-radius equilibrium below a threshold 
surface energy. 



HPq 1 sec 6> p . Hence, Fbuik/(^VX) favors small bundles, 
acting against lateral bundle growth. 

As shown in Fig. [5] , below a critical value of surface 
energy, E c , the minimal free energy state is characterized 
by a finite bundle size. The critical surface energy derives 
from the solution to 



/bulk(oo) — /bulk(-Rc) — R, 



dfbu 



ik 



dR 



L 



Rr 



(48) 



where /buik(-Rc) = F^ u \^/N and R c is finite-size of bun- 
dles at the point of equilibrium with the bulk assembly 
(R — > oo). For cohesive bundles this yields 



R el sin 2 9 V ; E^ las = S? las sin 6 9 V ( S? las = Bk RclPo 



:-52 



while for osmotic bundles, it gives 



(49) 



R° 



2R os ^J '(sec 9 P - l)/3 for R° s < R os sin 
(sec^-1) 3 / 2 (ST^). 



When R° s < R 

os sin 9 p the optimal bundle remains iso- 
metric at the critical surface energy, which holds only 
for 9 P < 0.71. For larger helix angles, the packing and 
size of osmotic bundles must be determined numerically. 
We may estimate the thermodynamic properties in this 
limit, however, by approximating equilibrium bundle size 
to be at the critical value, i?i so = R os sin9 p . Since 
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^tot (Ros sin 6) this overestimates the free energy of the 
finite-size bundle, this approximation gives a lower bound 
on the critical surface energy, 

£° sm > yi?o S sin^[/ bulk (oo)-/ bulk (i? os sin^)]. (51) 

In Figure [6] we plot the phase diagrams for models 
of cohesion- and pressure-driven condensation of heli- 
cal filaments, where a critical surface energy S c sepa- 
rates a phase of finite-size bundles from bulk assembly 
(in the isomorphic packing). From eqs. (|49|) and (jSTj) 
it is easy to show that for small helix angles that criti- 
cal surface energy increases as power law for both mod- 
els: E° las - 6* and £° sm - 0|. This gives the intu- 
itive result that as the preferred filament approaches a 
straight geometry, thermodynamic costs associated with 
packing frustration must vanish. Minimization of the 
surface energy (for any value of S) for bundles of rod- 
like filaments results in aggregates of unlimited lateral 
size, the result of standard arguments for 2D aggrega- 
tion. In the large helix angle limit, the costs of inter- 
filament packing of the isomorphic, R — ¥ oo limit increase 
in both models. Consequently, the critical surface energy 
needed to offset the bulk costs of isomorphic state also in- 
creases indicating a broader thermodynamic stability of 
finite-sized bundles. Hence for cohesive bundles we find 
~ Rf as fj> sin 4 6p ~ sin 6 # p , and a maximal stability 
for bundles at P — > 7r/2, or the limit k 3> r. We find 
a profoundly different result for the 9 p — > ir/2 limit of 
osmotically-driven bundles, as E° sm diverges due to di- 
vergent cost of isomorphic packing of the "ring-filament" 
limit. 



IV. DISCUSSION & CONCLUSION 

In this paper, we studied a continuum model for the 
formation of bundles of intrinsically-curved, helical fila- 
ments. A key aspect of this problem is a fundamental 
frustration between the curved shape of individual fil- 
aments and the tendency for even filament spacing, or 
maximal density, in the cross section. As implied by 
cq. ([1]), isometric packing requires filament distortions, 
whose costs increase rapidly with cross-section thickness. 
In comparison, the thermodynamic costs of imperfect 
packing in both models is relatively insensitive to bundle 
radius. While the (per filament) cost of filament distor- 
tion in isometric bundles grows as ~ k 4 R 2 , according to 
eqs. ([271 ) and the inter-filament packing cost in 
isomorphic (per filament) is independent of radius. This 
difference in size sensitivity underlies the generic evolu- 
tion from isometric packings for small bundles (R — > 0), 
where mechanical distortion is negligible, to isomorphic 
packings for large bundles (R — > oo), where filament dis- 
tortion becomes prohibitively expensive. 

Despite the common tendency towards constant shape 
packings in larger bungles, we find a markedly different 
sensitivity of inter-filament packing to lateral bundle ra- 



1.2 - 




FIG. 6: The phase diagrams of cohesively-driven and 
osmotically-driven assembly of helical filaments, respectively 
in (a) and (b), in terms of helical angle of filaments and sur- 
face energy bundles. Regions of self-limited (bundles) and 
bulk assembly are indicated. 



dius for cohesive bundles and osmotically-condensed bun- 
dles, as evidenced by the i?-dependence of optimal pack- 
ing angles shown in Figs. [3] and 01 In cohesive bundles, 
packing evolves smoothly with R for all helix angles and 
mechanical properties, exhibiting a non-isometric pack- 
ing (i.e. 6? e i 7^ 0) for any finite radius. In contrast, when 
the cost of filament torsion is negligible (either f3 — > 
or 0p — > 7r/2), in osmotically-condensed bundles, we find 
that bundles "lock-in" to the isometric packing over fi- 
nite range of small radii, R < R osra shi9p, exhibiting a 
continuous, second-order transition to the non-isometric 
packing beyond this critical size. The difference in in- 
dependence of the packing in the two models derive ul- 
timately from a different "stiffness" to interfilament dis- 
tortion from the isomorphic state, <5(|Aj_|) ~ —9 2 . In 
the elastic model this leads directly to a packing angle 
dependence, F c i as ~ |d(|Aj_|) ~ 9 4 , indicating a "soft" 
^-dependence of packing cost that is generically weaker 
than the ^-dependence of -Emoch driving bundles toward 
the 9 — > 9 P isomorphic state. This can be contrasted 
with the "stiffer" 0-dependance of packing free energy in 
osmotic bundles, F c \ &s oc AV ~ <5(|Aj_|)| ~ 9 2 , a com- 
parable dependence to the bending costs, which implies 
a finite size range over which the isometric packing is 
stable. 

The evolution from isometric to isomorphic packing 
has important consequences for the equilibrium size of 
self-assembled bundles of helical filaments. In ref. [||, it 
was shown that constraining bundles of helical filament 
to the isometric packing leads to self-limited lateral radii 
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of equilibrium bundles. Hence, this system belongs to 
an unusual class of associating molecular systems along 
with self-twisted bundles of chiral filaments 0, H, I2H - 
[27j . where the equilibrium dimensions of the condensed 
state are finite, but mesoscopic — much larger that fila- 
ment diameter — despite the finite range of the inter- 
actions. Importantly, the present study demonstrates 
that even when inter-filament packing is free to relax 
from the isometric packing, the self-limiting behavior of 
helical filament bundles is retained over a broad range 
of cohesive conditions. From Fig. [5J we note again a 
significant difference between the stable range of self- 
limited bundles in cohesive and osmotic bundles, char- 
acterized the large-# p limit of the E c , the surface en- 
ergy where finite-size bundles are in equilibrium with 
bulk (infinite- R) assembly. While finite-diameter cohe- 
sive bundles reach a maximal range of stability in ring- 
filament limit, Ef(6 p ->• tt) = E cl , E° sm diverges in this 
limit, indicating that the stable range of finite bundles 
becomes infinite in the k/t»1 limit, where the costs of 
packing frustration are maximal. 

We conclude with a brief discussion of relevant phys- 
ical parameters for a model system of helical filaments, 
the bacterial flagella. Owing the precise control of heli- 
cal geometry offered through the multiple polymorphic 
shapes of flagella HH , the condensed states of recon- 
stituted flagella provide an ideal system for exploring 
the thermodynamics of curved-filament assembly. Re- 
cent experiments have exploited the intrinsic shape of 
reconstituted flagella in colloidal suspensions as a means 
to study study the unusual liquid crystalline mesophases 
of spiral-shaped filaments [20|. Here, we consider the 
properties of bundles of flagella that may be formed in 
highly-concentrated solutions based on our analysis. For 
simplicity, we focus on the case of bundles of flagella 
osmotically-condensed in the presence of a osmolytes suf- 
ficiently large to be excluded from the interstitial regions 
of the dense bundle core. Of particular importance is 
the size scale, R os in eq. (fl3")l . which determines both 
the range over which bundles are isometrically packed 
(R < R os sm9 p for j3 = 0) as well as size range of finite- 
size bundles in equilibrium according to eq. (l5"Tj) . We re- 
cast the definition of this length as R os = diJl/Hmcch) 1 ^ 2 , 
where 



n 



mcch 



Bk, 4 



(52) 



constitutes a pressure scale determined by the deforma- 
tion cost of filaments (we have dropped numerical fac- 
tors of order unity). Force-extension measurements of 
flagella under flow suggest a bending modulus, B = 
3.5 pN//zm 2 (28|. Using the preferred curvature of the 
normal and coiled polymorphs, K norma i = 1.3 /im -1 and 
Kcoiied = 1-8 jUm , we find predictions for the char- 
acteristic osmotic pressures, Tl(normal) m 10 Pa and 
H(coiled) « 38 Pa. Compared to the typical scale of 
osmotic pressure — of order atm ss 10 5 Pa for few per- 
cent of 8000 MW PEG — the modest pressure scale set 
by flagellar bending cost suggests that the characteristic 



bundle size is mesoscopic, orders of magnitude larger fila- 
ment diameter, i? osm ~ 10 2 d. Since sin6* p s» 1 for helical 
angles of these flagellar polymorphs — 9 p (normal) ~ 31° 
and 9 p (coiled) ~ 77° — our theory suggests that pack- 
ings remain isomorphic for bundle radii up to the range 
of microns. 

Turning now to the range of self-limited lateral assem- 
bly, we consider the surface energy of bundle under the 
influence of osmotic pressure, E = Tlr, where r is the 
effective size of that "depletion zone" (roughly, the os- 
molyte radius) that excludes osmolytes near the bundle 
surface [30[ ■ From the critical surface energy of eq. (jSTj) , 
we expect finite diameter bundles above a critical pres- 
sure, n c w (r/d) 2 n moch (sec6» p - 1)~ 3 , again suggesting 
the threshold osmotic pressures are nominal, in the range 
the 1 — 10 Pa. Thus, we expect reasonable experimen- 
tal conditions to be in the range of self-limited bundle 
sizes, for which the balance between surface energy and 
filament distortion in isometric bundles determines the 
equilibrium dependence of bundle size on osmotic pres- 



R/d; 



(n/n moch ) 1/3 , 



(53) 



which is of order 10 for n « 1 atm. These estimates 
suggest a strong (measurable) pressure dependence self- 
limited diameter of flagellar bundles in the range of weak 
osmotic pressures n > n mec h. 
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Appendix A: Derivation of Shape Variation 

Here, we derive the shape variation of filaments in the 
family of bundles described by eqs. (U])-©. The unit 
tangent, T(x), of filaments at x follows straightforwardly 
from eq. (|4|) 



T(x) 
where 



ds (x) 



ds 



(T + ST) ~ T(l - T • ST) + ST, (Al) 



ST = <9 s x = (ft + ft )p x x- ft p (N- x)p, (A2) 



and 



(ST ■ T) = [(ft + ft ) sin 9 - cos 9fl p ] (N • x). (A3) 
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ds(x) 
ds 



ST\ ~ 1 + T-5T, 



(A4) 



reflects mapping of arc-length elements at x to the central 
filament arc-length. Recall in this notation objects that 
appear as explicit functions of x denote geometry quan- 
tities at position x in the packing plane, and objects not 
written as functions of x refer to geometric quantities re- 
lated to the central filament at x = 0. Note also that 
we are working in the limit of narrow bundles where it 
is sufficient to consider first-order in |x| correction to the 
geometry (relative to the central filament). 

The filament curvature derives from the magnitude of 



(A5) 
(A6) 



relative to the central filament. The normal correction 
follows from, 



5N 



N 



1 dT(x) 
k(x) g?s(x) 
k _1 (B • d s 5T)~B — T(N • 6T), (All) 



riT(x) _ kN(1 - T ■ ST) - Td s (T ■ ST) + d s ST 
ds(x) ~ 1 + T • ST ' 

from which we derive 

k(x) - k ~ N- d s ST ~ 2k(T ■ ST), 

where 

8 s st ~ -(ft + n fx + n p (n + n )W x x 

+ ft p ft(v' ■ x)p - ft 2 (u' ■ x)u'. (A7) 
Using the fact that u' = N is perpendicular to p we find, 
N -d s ST = -[(ft + ft ) 2 + ft 2 ](N -x), (A8) 
which when substituted into eq. (|A6[) yields eq. (??). 

6k = -«T X [fft+ft ) 2 +ft p ] -2 [(ft+ft o )sin0-cos0ft p ], 

(A9) 

where k(x) — k = Sk(kN ■ x). 

The variation of torsion follows from 

r(x) = -N(x) • — t(1 — T ■ 6T) — N • d s <5B ± , 

(A10) 

where 6N± ~ N(x)-N and <5Bj_ ~ B(x)-B are the first 
order variation to normal and bi-normal, respectively, 



where we used N(x) • T(x) — 0(|x| 2 ) in the second line. 
Similarly, defining ST±_ = 5T - T(T • ST), we have the 
correction to the bi-normal, 

5B±_ ~ T x SN_l + 5T ± x N 

= -N(B • dsST)^ 1 - T(B • ST). (A12) 



The N component of the derivative of eq. (IA12[) can be 
written as 



N • d s SB ± = -(N • a s 5T) - (B • d^ST)^ 1 - k(B ■ ST). 

(A13) 

From eq. (|A3[) we have 



(B-5T) = [cos0(fi + Oo)+sin6l^p](N-x), (A14) 
and differentiating eq. (| A7|) we find 

(B ■ d 2 ST) = | - (Q + fl ) 2 [ cos 0(fi + n ) + sin 0fi p ] 

- rj p fi(fi + Q ) sin0 - sin0^ p ^ 2 - rj 2 r|(N • x). 

(A15) 

Substituting these into eq. (IA10[) we find 



n+nn) 2 +nl~\ 



St = [(O+fio) sin 0— cos 0f2 p ] — k 2 |— r [(si+s2o 

+ [k 2 + (ft + fi ) 2 ] [ cos 0(fi + n ) + sin 0ft p ] 
- ft p [ sin 0(ft + ft )ft + sin 0ft 2 + ft p r] j, (A16) 



where r(x) — r = <5t(kN • x) 
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